Exact statistical properties of the zeros of 
complex random polynomials 



P.J. Forrester and G. Honner 
Department of Mathematics and Statistics, 
University of Melbourne, Parkville, Victoria 3052, Australia 



The zeros of complex Gaussian random polynomials, with coefficients such that 
the density in the underlying complex space is uniform, are known to have the 
same statistical properties as the zeros of the coherent state representation of one- 
dimensional chaotic quantum systems. We extend the interpretation of these poly- 
nomials by showing that they also arise as the wave function for a quantum particle 
in a magnetic field constructed from a random superposition of states in the low- 
est Landau level. A study of the statistical properties of the zeros is undertaken 
using exact formulas for the one and two point distribution functions. Attention 
is focussed on the moments of the two-point correlation in the bulk, the variance 
of a linear statistic, and the asymptotic form of the two-point correlation at the 
boundary. A comparison is made with the same quantities for the eigenvalues of 
complex Gaussian random matrices. 



1 Introduction 

The zeros of random polynomials have received recent attention as providing the universality 
class for the zeros of the coherent state representation of one-dimensional chaotic quantum 
systems Moreover, denoting the coefficients by aj so the polynomial reads 

N 

P{z) = ^ajZ^, (1.1) 
i=o 

it has recently been realized that in the particular case that the real and imaginary parts of aj 
are independent Gaussian random variables with mean zero and standard deviation aj, which 
we write as 

= N[0,crj] -hiN[0,crj], (1.2) 



the zero distribution is a solvable model |0 - |lj]. Thus the probability density function (p.d.f.) 
can be computed in terms of the aj 0, as can the n-particle distribution functions for each 
n = 1,2,... m. 
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(1.3) 



Our study is motivated by the explicit form of the zeros p.d.f., 
where 

Cj := ^ ^pi ' ' ' ^Pj (1-4) 

l<pi<---<Pj<N 

and Co := 1- As pointed out in Q, the factor 

n \zk-Zj\^ = e-^^^<^'°^^''-''^ (1.5) 

^<j<k<N 

can be interpreted as the Boltzmann factor for a two-dimensional classical particle system 
with pairwise logarithmic repulsion at inverse temperature f3 = 2. The particles are pre- 
vented from repelling to infinity by an extensive many body interaction with Boltzmann factor 
{J2jLo l^iP/^AT-j)"^- Furthermore, with the particular standard deviation 

(1.6) 



(j!)l/2 



the density of zeros is, to leading order, uniform inside a disk of radius vA^ and zero outside 
this disk Q (see also Section 5). This latter feature is shared by the eigenvalue p.d.f. for 
complex random matrices with each element independent and chosen from the complex Gaussian 
distribution N[0, l/\/2] + iN[0, l/\/2], which is proportional to 



ne-s.-i^^i n \^k-z,\'. (1.7) 

j=l l<j<k<N 



Let us consider further the eigenvalue p.d.f. ( |1.7[ ). Some insight into the behaviour of the 
density can be obtained by noting that for \z\ < \/]V, (^]^) is proportional to the Boltzmann 
factor of the two-dimensional one-component plasma (2dOCP) at the coupling f3 = 2. The 
factor exp(— J2jLi kjP) there results from the potential energy between a uniform background 
of charge density — l/vr and A^ unit charges interacting via the logarithmic potential, so the fact 
that the density has support in a disk of radius ^/N can be interpreted as the system having a 
preference for charge neutrality. Moreover, the plasma system is an example of a two-dimensional 
Coulomb system in its conductive phase, and so exhibits universal properties for the behaviour 
of the two-point correlation, both in the bulk and at the boundary \z\ ~ \/]V (see e.g. 16|). 



These special forms in turn imply a special form for the variance of a linear statistic |17]. The 



eigenvalue distribution of course must also exhibit these Coulomb properties. Because the zeros 
p.d.f. has the factor (|1.5| ) in common with (|1.7| ), and with the choice ( |1.6D has the density in 
common too, we are led to investigate for the zeros distribution the two-point correlation in the 
bulk and at the boundary \z\ ~ ^/N, and the variance of a linear statistic, to see if they exhibit 
universal Coulomb like properties. 

The plan of the paper is as follows. In Section 2 an interpretation of the complex Gaussian 
random polynomial with variance (|1.6|) will be given as a random wave function associated with 
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a superposition of lowest Landau states in the plane. A similar interpretation will be given to 
the complex Gaussian random polynomial with variance of the coefficients 

considered originally in Q, except this time the Landau states are those for a particle confined 
to the surface of the sphere with a perpendicular magnetic field. In Section 3 the known |^ 
general formulas for the one and two-point distributions are revised, and these are specialized to 
the cases and (|l.8|). The investigation proper is carried out in Sections 4 and 5. In Section 
4 we study bulk properties of the two-point correlation function, and show how the formulas 
obtained for its moments have consequence with regard to fluctuation formulas for the number 
of particles in a region A as |A| — > 0, and the variance of a linear statistic as a function of in 
the sphere case ( |l.8D . The behaviour of the density and two-point correlation at the boundary of 
support in the disk case ( |1.6D is the subject of Section 5. The paper ends with a brief summary 
of our results. 



2 Quantum mechanical interpretation 

A number of studies Q-Q have related complex polynomials to the wave function of chaotic 



quantum systems. More precisely, complex random polynomials with variances and (1.8) 
arise as the Bargmann or coherent state representation of one-dimensional chaotic wave functions 
and chaotic spin states respectively. In this section we point out that these complex random 
polynomials also occur as random superpositions of lowest Landau level states for a charged 
quantum particle, in a plane and on the surface of a sphere, in the presence of a perpendicular 
magnetic field. 

Consider first the planar problem. Let B denote the strength of the magnetic field, and 
suppose the vector potential A is chosen in the symmetric gauge so that 

Then it is well known (see e.g. fl^ ) that the lowest Landau level (i.e. ground state) consists of 
any wave function of the form 

e-l^l^/^^V(^-), (2.1) 



where z := x + iy, I := yJfic/eB (e denotes the charge of the particles) is the magnetic length 
and / is analytic in z. Furthermore the states ( |2.1D can be separated by requiring that they 
be simultaneous eigenstates of the quantum mechanical operator corresponding to the square of 
the centre of the cyclotron orbit. This operator has eigenvalues -R^ = (2n -|- 1)1'^ (n = 0, 1, . . .) 
with corresponding normalized eigenfunctions of the form (|2.1|) given by 



3 



(here the units have been chosen so that 2P = 1). A state 0(r) consisting of the first + 1 of 
these states has the form 



^('^ = tA72^"''''^'p(^")' p(^) ■■= E (2-3) 
A/vrV^ Vn! 

with J\f := {J2jLo IcKnP)"^''^- Choosing the state to be random by selecting each coefficient at 
random from the complex Gaussian distribution N[0, 1] +iN[0, 1], we see immediately that p{w) 
is a complex Gaussian random polynomial with variance ( |1.6D . 

For the sphere problem, the magnetic field is due to a monopole, which must satisy [^0| the 
quantization condition B = N'hc/2eE? , N = 1, 2, . . .. Furthermore the Hamiltonian operator H 



can be written |f9| 



where the components of L obey the canonical commutation relations for angular momentum 



[Lx,Ly] — ihLz, [Ly,Lz] — ifiLx, [Lz,Lx] — ifiL 



y 



Hence the allowed values of are l{l + I = 0, ^, 1, |, . . . and so for ( p.4[ ) to be positive 

definite the smallest allowed value of I is N/2, which thus corresponds to the lowest Landau 
level. The states in this level can be distinguished by seeking simultaneous eigenfunctions of 
and H. This shows that there is a + 1 fold degeneracy, which is spanned by the orthogonal 
functions 

^-(^'■^^ = {^lS^(Ar/+m)}'^'"''^'^'""''^'"'"'"^^^^^ = -Ar/2,-Ar/2 + l,... ,Ar/2, 
where 

2 2 

denote the Cayley-Klein parameters and (9, (p) are the usual spherical coordinates. Forming a 
linear combination of such states gives the state 



1 ^Ar/2 + 1^1/2^ -isin0/2)^e-^*/Ve^<^cot-), p{z) = V (^Y^nZ^- 



M\ 47ri?2 



If this state is chosen at random by specifying the an as belonging to the random complex 
Gaussian distribution N[0, 1] + iN[0, 1] as in ( |2.3| ), then the polynomial p{z) is a complex random 
polynomial with variance (ITsI) . Further more, the mapping 



z = e'^ cot - (2.5) 

represents a stereographic projection from the north pole of the unit sphere to a plane passing 
throught its equator, showing that p as a function the spherical coordinates {9,(j)) has A^ zeros 



on the surface of the unit sphere. Such a factorization of a spin state is due to Majorana [21|. 
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3 The distribution functions 



As mentioned in the Introduction, the general n-point distribution function for the zeros 
of the random polynomial (|1.1|) can be computed exactly in the case that the coefficients are 



complex Gaussian random variables drawn from the distribution (1^) Q. Our interest is in this 
explicit form in the cases n = 1 and n = 2. 

The case n = 1 corresponds to the density of zeros. This one point distribution function has 
the explicit form 1^, 14 1 



log{p{z)p{z)) , 



(3.6) 



dzdz 

where the averages are with respect to the distribution of the coefficients. Here r := (x, y) with 
X + iy := z. Performing this average in the case that the distribution of the coefficients aj is 
given by ( L£ ) gives 
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dzdz 



loge(|z|2;iV) 



where 



3(n;n) := ^ ^. 



p=0 



Notice that with z fixed e(|2;p; N) — > e'^'^ as — > 00, so the formula ( ^.7] ) gives that 



P{i){r) = - 



(3.7) 
(3.8) 

(3.9) 



in the bulk. 



Repeating the calculation which led to (^]^) with the coefficients given by (1.8) shows that 
the corresponding density is 

92 



7r/>(i)(0 



dzdz 



log 1 + \z\ 



N 



(3.10) 



Due to the physical origin of the variances ([1.6|), it is natural to perform a stereographic projec- 
tion onto the sphere according to the mapping (|2.5D so that 



(3.11) 



P(i){r)dxdy = P{i)i^) ^^^^\l\2y 



where p{s) refers to the density on the unit sphere with surface element dS at the point s 
determined by (2^) and 4/(1 + \z\'^)'^ is the Jacobian for the change of variables. This gives 

& 113 

N 

Pii)is) = j-, (3.12) 



thus showing that the density is uniform. 

Not surprisingly, the complexity of the explicit form of the n-point distribution increases 
with n. We will not go beyond the n = 2 case, for which the distribution has the explicit form 

per(C - B^A-^B) 



P(2)(ri,r2) 



detA 



(3.13) 
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The notation per denotes the permanent, which for a 2 x 2 matrix is defined by 



per 



and with pj := p{zj), p'j := p'{zj), (j 
A, B and C are defined by 



" ^ ■=ad + hc, (3.14) 
c d \ 

1, 2 and the dash denotes differentiation, the matrices 



A := 
Writing 



(PiPi) (P1P2) 
{P2P1) {P2P2) 



B 



A-' 



detA 



{Pip'i) {PiP'2) 

{P2p'i) {P2P2) 

«22 —012 

-021 «11 



C :-- 



1 



(p'lP'i) {p'iP'2) 
{P2P1) {P2P2) 



(3.15) 



detA 



-A 



D 



we see that ( 3.13 ) can be written in the equivalent form 



TT' 



'p{2){ri,f2) 



per(dC - B^A^B) 
d3 ^ 



d = detA. 



(3.16) 
These are 

(3.17) 

where e{u;n) is defined by ( |3.8| ), which in turn fohow from the definitions. At this stage we 
will not attempt to further develop ( |3.16| ); this task will be undertaken in Section 5 when the 
^ 00 limit in the neighbourhood of the boundary \z\ = y/N is computed. However we will 
note the previously computed ||9| A ^ cxd limiting value with ri and f2 fixed. This limit gives 
the 2-point distribution function in the bulk, and has the explicit evaluation]^ 



Consider now the evaluation of the matrix elements in ( ^.15 ) in the case of ([L 
made explicit by referring to the formulas 

{PiPj) = e{ziZj;N) {p'-pj) = Zje{ziZj;N - 1) 

{PiPj) = Zie{ziZj] A - 1) (p-p^-) = ZiZje{ziZj;N - 2) + e{ziZj]N - 1), 



P(2){n,r2) 



1 



/i(|fi-f2|2/2), Kx):-- 



(sinh^ X + x^) cosh x — 2x sinh x 



(3.18) 



sinh X 

In the case of the variances (O) we also require knowledge of the matrix elements in ( |3.15| ). 
A simple calculation gives [0] 



{PiPj) 



{l + zfz.f 



{p'iPj) = Nzj{l + ZiZj) 



N-l 



{piP'j) = Nzi{l + z,Zjf~^ (pWj) = N{1 + NziZj){l + z,zj)^~\ (3.19) 

In Section 4 we will have application for the corresponding evaluation of ( |3.16| ) in the finite 
system, when projected onto the sphere according to the transformation ( |2.5| ). By substituting 
( 3.19| ) into ( 3.15| ), and then the result into ( |3.16| ), expanding the permanent according to ( p. 14 ) 
and performing some elementary but tedious manipulations we find 



P(^)(S1,S2) 



P(2)(^l,S2) - P(1)(S1)P(1)(S2) 
Ar2 fN-4 



167r2 1 - / 



27V 



(1./ 



2\2, 



1-/ 



2N 



2N\2 



(3.20) 



^ In Q the factor of 1/2 in the argument of /i(|ri — r2p/2) is missing; this has been subsequently 
corrected in dj. 
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where 



(3.21) 



4 The distribution functions in the bulk and fluctu- 
ation formulas 

4.1 The two-point function in the bulk 

In this section properties of the bulk two-point distribution function for the 2dOCP at /3 = 2, or 
equivalently for the eigenvalues of complex random matrices specified in the Introduction, will 
be recalled and contrasted with the corresponding properties of the bulk two-point distribution 
function ( p. 18 ). Now, for a bulk density P(i)(r) = 1/vr, the bulk two-point distribution function 



for the 2dOCP at /3 = 2 has the exphcit form |2|, |23| 



P0CP(2)(ri,r2) = ^(l - e-l^^i-^^^l'), (4.1) 



or equivalently the truncated distribution has the explicit form 

i1 (t, Tr,^ = — Lp-l'^i-^al 



pScP(2)(n,r2) = --e-l^^-^^^l'. (4.2) 



vr 



One feature of (|4.l| ) is that it is a strictly negative monotonic function and is its own large 
separation asymptotic form, 

Pocp(2)(^,0) ~ --^e II. 

In contrast ( p.l8| ) reaches a maximum before approaching asymptotically approaching zero from 
above (see e.g. the plot in Q). The large separation asymptotic form is computed from ( p. 18 ) 
to be 



Pf2)(^"',0), ~ |rlV 



The zeroth and second moments of /Oqcp (2) exhibit particular universal properties of Coulomb 
systems. The zeroth moment is evaluated as 

1 



which can be written in the equivalent form 

/ , (pScp (2) i^^ 0) + P{1) (f^^i^) df = 0. (4.3) 

In this guise the zeroth moment formula can be interpreted as an example of a charge sum rule 
for Coulomb systems (see e.g. p5[]) which says the total charge of a fixed charge (here a charge 
at the origin) and its screening cloud (represented by /OocP(2)(^' ^)) zero. 

The second moment of ( [4.1[ ) is evaluated as 

[ M'plcPi2)ir,0)dr = -l. (4.4) 
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By invoking a linear response relation (see e.g. |15, |l^) this can be shown to be equivalent to the 
statement that the system perfectly screens an external charge density in the long wavelength 
limit. Thus suppose the 2dOCP for general (3 is perturbed by an external charge density ee^^'"^ . 
The potential energy of the system is then changed by an amount 

5U := — e / dr' / dr \og\r — r'\n(i\{f')e^^'^ 

(4.5) 



\k\ 



where represents the microscopic density at the point r, ^i(i)(A;) its two-dimensional 

Fourier transform, and to obtain the second equality the result 

/■ ^ ^ 27r 
— / dr loglfl = -^r- (4.6) 

from the theory of generalized functions has been used. Suppose now that the microscopic 
charge density at the point r, which for the OCP is just n(i)(r), is observed. The linear response 
relation gives that 

(n(i)(f)), - (n(i)(f))o = -/3(n(i)(r)<^C7)^, (4.7) 

where the subscript (e) indicates the system in the presence of the perturbation. Now, a charac- 
teristic of a Coulomb system in its conductive phase is that it will perfectly screen an external 
charge density in the long wavelength limit. This means that 

(n(i)(r)),-(n(i)(f))o ~ -ee*^: 

Substituting this in (^]^), and using the translational invariance of the system in the bulk, we 
see that 

S{k) - S{r) := pI^{t,^) + p5{r) (4.8) 



This implies the charge sum formula ([4.3| ) as well as the second moment (Stillinger-Lovett) sum 
rule 

df|rf p5cp(2)(^"'>0) = (4.9) 
which reduces to (O) at (3 = 2. 



Let us now compute the zeroth and second moments of ( 3.18| ), and compare their values 
with those for ( [4.1| ). Now, whereas the computation of the moments in the case of ( [4.1D is 
elementary, it is not so simple to do likewise for ( 3.18| ). Fortunately the function h{x) in ( p. 18 ) 
can be written as the second derivative of an elementary function. Thus by explicit calculation 
we can check that 

1 d'^ 

h(x) = — — rx^cothx, (4.10) 

2 ax^ 

or equivalently 

/i(x) - 1 = i-^(x^(cothx - 1)). (4.11) 
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Making use of ( 4.11| ) we find that again the charge sum rule (O) is satisfied, whereas for the 
second moment we find 

J^^dr\r\^pl^{f,6)=0. (4.12) 

Although we have know that when viewed as a classical particle system the p.d.f. ( |1.3D for 
the zeros distribution contains many body forces in addition to the pair potential term (|1.5| ), 
it is of some interest to compare the results ( |4.3D and ( [4.12 ) to those for a two-dimensional 
classical fluid with purely pairwise interactions. For long-range potentials with Fourier transform 
v{k) ~ c|/c|~", a > 0, we can generalize the perfect screening argument which led to ([4.^), this 
time obtaining 

S{k) ~ ^. (4.13) 

\k\~>Q CfJ 

Setting k = gives (^^) in all cases. However, as is well known (see e.g. [^), it is only for 
even values of a that (|4.13| ) is compatible with a fast decay of S'(r). In the Coulomb case a = 2, 
( 4.13| ) then gives the second moment condition ( [4.91 ). The next smallest even value is a = 4. 
Then (4.13) indeed implies that the second moment vanishes, just as observed in ( 4.12| ) for the 
complex zeros. However a pair potential in two-dimensions for which the Fourier transform 
behaves as \k\'^ for \k\ 0, has the large \r\ behaviour in real space proportional to —\r\'^ log |r|. 
It is unlikely that a classical system with such strong repulsion can be thermodynamically stable, 
even with a neutralizing background. 



4.2 Fluctuation formulas in the infinite system 

In general the properties of the two-point function of a many particle system, in particular the 
value of its moments, have consequences regarding fluctuation formulas for linear statistics. We 
recall that a linear statistic is any observable (function) of the form A = X]jLio(^)- It is 
easy to see from the definitions that the two-point function S{f) ([4.8[) is directly related to the 
fluctuation (variance) of a linear statistic via the formula 

Var(^) := {A^) - {Af 

dfia{fi) / df2 a{f2) S (r) ■ (4-14) 



The linear statistic specified by a(r) = x^i^ where XA denotes the indicator function of 
the domain A: xa(^) = 1 if r G A, XA(j^) = otherwise, measures the number of particles in 
the region A. In this case Var(yl) is the number variance, to be denoted S2(A) which specified 



the fluctuation from the mean of the number of particles in A. It is a well known fact [26| that 
whenever A can be generated by a dilation from a flxed domain Aq, and the charge sum rule 
(4.3) holds, the number variance has the asymptotic behaviour 

hm = -- f , |r>5)(f,0), (4.15) 

|A|-^oo \oA\ TT J-Ri > 

assuming the integral exists. Here, with A a two-dimensional domain, \dh\ denotes the length 
of the perimeter of A. Thus the fluctuations are strongly suppressed in comparison to a com- 
pressible fluid for which S2(A) is proportional to |A|. 
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In the case of the complex zeros the two-point function is specified by ( 3.18 ). Making use of 
the formula ( |4.11| ) the first moment in (4.15) can be computed exactly. For this purpose, after 
introducing polar coordinates, we integrate by parts twice and make use of the definite integral 



1 



^C(-) 



Substituting the result in (|4.15D gives 

S2(A) 

lim 



87r3/2''4^ 



0.0586436 • 



(4.16) 



|A|-oo |M| 

where (^{s) denotes the Riemann zeta function. In the special case that A is a disk of radius 
r, Prosen ||l^ has obtained a one dimensional integral representation of S2(A), and from this 
obtained a numerical evaluation consistent with ( 4.16| ). 

Another class of linear statistics of particular interest with regard to the eigenvalues of 
random matrices (see e.g. |17] and references therein) is when o(r) varies on macroscopic length 
scales relative to the spacing between particles. This can be achieved by first supposing o(r) 
varies on the length scale of the spacing between zeros, and then making the replacement a(r) i— > 
a{f/a) where a ^ 1. Now, writing S{r) in ( 4.14 ) in terms of its Fourier transform we have the 
formula 



Var(yl) 



1 



Making the replacement a{r) 



(27r)2 jR2 
a{r/a) then gives 

2 



dr\~a{k)\^S{k). 



Var(yl) 



a 



(27r)2 7r2 



dr\~a{k)YS{k/a). 



(4.17) 



Since a ^ 1 we see that the small |fe| behaviour of S{k) determines Var(^). For the OCP 
S{k) exhibits the quadratic behaviour (|4.8D and so Var(yl) is formally independent of a in the 



macroscopic limit [16|. For the complex zeros, as already noted the fact that the zeroth and 
second moments of S{r) vanishes implies that 

ck\ 



S{k) 



(4.18) 



Furthermore, from the definitions it is easy to see that c is related to the fourth moment of pj^-^ 
by 



1 



- I^Jfr P(2)(f,0). 



R 



Using ( [4.11| ) this integral can be computed exactly, giving 

1 



87r 



C(3) 



(in fact all coefficients in the power series expansion of S{k) have been calculated in 

Substituting ( [4.18| ) in ( |4.17D we see that 

/ c 



Var(A) 



V(2vr)2 



R2 



dk\d{k)\'\k\' 



assuming the integral exists. (Note that the integral does not exist for a(r) 
with the distinct behaviours of ( |4.15D and ( 4. IS] ).) 



(4.19) 
XAo, in keeping 
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4.3 Fluctuation formulas on the sphere 



The fluctuation formula ( [4.19|) applies to the infinite system. Also of interest is the situation 
of a fixed volume, and thus the density being proportional to as the number of particles 
increases. This is the natural setting for the zeros of the complex random polynomial with 
variances specified by ( |1.6D projected onto the unit sphere. We know from ( |3.12 ) that the zeros 
are uniformly distributed. The question of interest is the behaviour of Var(yl) as a function of 

In fact this question, and its generalization to variance of linear statistics for zeros of mul- 
tidimensional polynomials projected on to the surface of higher dimensional spheres for which 



the zero density is uniform, has been addressed recently by Shiffman and Zelditch |14]. They 
proved that for a statistic with a{6, (j)) smooth on the unit sphere. 



Var(A) = 0(1) 



(4.20) 



as A — > oo. This bound was also established in the higher dimensional cases. The question 
now is if the 0(1) behaviour is asymptotically exact. In light of the analysis for the infinite 
system presented in the previous subsection, we would expect that asymptotically exact 0(1) 
behaviour corresponds to S{k) ~ ci|A;p in the infinite system. This can be seen from ( 4.17| ) with 



a 



identified with A. But S{k) for the complex zeros has the small \k\ behaviour (4.18) which 



is quartic in Identifying with A in ( 4.191) then suggests that 



Var(^) ~ -. 



(4.21) 



We have not been able to prove this assertion in general, but it can be demonstrated in a specific 
example, using a combination of analytic and numerical analysis. 



According to the formula ( 4.14 ), our task is to compute 

Var(^) = I dSi / a(fi)a(f2)pL(fi, 7=2) + / dS a'^{r)p(^i){f), 



(4.22) 

m Jn ' ' Jn 

where dS is the differential surface element for the unit sphere, /0(i)(r) is given by ( 3.12| ) and 



P(2)i^i,r2) =pf2)(^i -^2) 



(4.23) 



is given by (|3.20[) . To proceed further, we suppose that in terms of the usual spherical coordinates 

a(f) = a(cos 9) = a{r ■ z) 
and thus depends only on the azimuthal angle. We want to use this functional dependence. 



together with the functional dependence displayed by (3.21) to simplify the double integral over 



the unit sphere in ( 4.22 ). 
First consider 



'1 := 



dSi a(fi • z)pj2){ri ■ h)- 



(4.24) 
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In this integral we change variables {9i,(j)i) {0,4') by rotating the sphere so that r2 ^ -2, 
where now 6 is measured as if the direction f2 were the z-axis and (f> = (p2- Then, in the new 
coordinate frame 

h = (0,0,1), 

fi = (cos sin ^, sin sin ^, cos 0), 



Thus 



= (— cos I 



!)2 sin 02, — sin (/)2 sin 02, cos 02) 



and so 
h = 



r ■ z = cos cos U2 — cos cp cos (p2 sm U sm U2 — sm (p sm 02 sm U sm ^2 

fi • f2 = cos 

j dS a(cos cos ^2 — cos (j) cos 02 sin sin 02 — sin </< sin (j)2 sin sin 02)P(2) (^^^ 0) . (4.25) 



If we now further specialize and choose 

a{fi ■ z) = fi ■ z 

then, after recalling that 

dS = sin 



(4.26) 



(4.27) 



we see that the integral over cj) in ( 4.25| ) can be carried out, leaving us with the single integral 



/i = 27rcos02 j d0 sm0 cos0p(^2)i^°^O). 
Substituting this simplification in 1 ^^ , together with (g]26|), (|]2^ and (|3]2|) we thus obtain 



TV 

Var(^fj • 



An J ^6*2 sin 6*2 cos 02 J (i0 sin cos 0/9(2) (cos 6*) 

N 
2" 



H / d0 sin cos 

2 Jo 



— — / dO sin cos 9 pjn) (cos 0) + — . 
3 JO 3 



(4.28) 



Although we have not been able to evaluate this integral analytically, the exact formula ( p.2[)| ) 
makes numerical integration straightforward. The results are given in Table 1, and are consistent 
with the expected behaviour ( [4.21| ). 



N 


Var(A) 


N 


Var(A) 


10 


0.1249 


80 


0.0193 


20 


0.0704 


90 


0.0172 


30 


0.0489 


100 


0.0156 


40 


0.0375 


110 


0.0142 


50 


0.0303 


120 


0.0130 


60 


0.0255 


130 


0.0120 


70 


0.0220 


140 


0.0112 
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Table 1. Numerical values of ( [4.28| ), truncated at the fourth decimal, exhibiting the expected 
0{1/N) decay of Var(A). 



5 Correlations near the boundary of support 

As remarked in the Introduction, to leading order the density of the complex zeros in the case of 
(|l.6|) is uniform in a disk of radius VN and zero outside that radius. The density of eigenvalues 
implied by the p.d.f. ( |l.7D have the same feature, and therefore so does the OCP at /3 = 2 
with the Boltzmann factor (1.7) taken to be valid throughout the plane (the harmonic potential 



between background and log-potential charge is only valid for the charge inside the background; 
outside the potential must be C + A^log|r| by Newton's theorem). As the two-point function 
of the OCP has a universal form at the boundary, it is of some interest to compute the form of 
the one and two point functions for the complex zeros, and compare their behaviour to that of 
the OCP. 



5.1 The density 



In the case that the variance of the coefficients of the random polynomial (1.1) is given by ( |1.6| ), 
we have noted that the density of zeros is given by ( |3.7| ). Here we want to analyze P[i){r) with 
r measured from the boundary of support of the density \r\ = \fN . This can be achieved by 
writing 

z = x + i{y-\^), (5.29) 

which sets the origin as the southern boundary of the support. According to (|3.7D the asymptotic 
expansion of e{\z\'^]N) with z a function of as specified by ( p. 29 ) is required. For this purpose 
we note from the definition ( |3.8| ) the easily verified formula 

,ni + N;z) 



e{z;N) 



T{N;a) 



T{l + N) 
and the asymptotic expansion [^] 

= - 1 — eri(n) + 0(—;=) 

r(i + A^) 2V ^ ^ 

Choosing u = (z — N) / V^N we see that 

e{z;N) = ie^(l-erf(^^) +0(^)), 

valid for (z — N)/\/2N bounded, which in turn implies 

e(|x + i{y - \/iV)|2; N) ~ ^e'^l' (l + erf (V2y) 



dt, 



(5.30) 



(5.31) 



(5.32) 



where z is given by ( 5.29| ). Hence, with the limiting form of P(i){r) denoted by P(i)(y), we see 
by substituting (|5.32| ) in (|3?7|) that 

1 5^ 

vr/3(i) (y) = 1 + 4 5^ l°g (l + erf (^^y)) • (5.33) 
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A plot of vr/9(i)(y), which is most expediently carried out by first explicitly calculating the 
second derivative to give 



, ;f2^e-V(l + erf(V2y)) + fe- 
^'^^^^'^ = ' - ^ (l + erf(V2.))^ ' ^'-''^ 

indicates rapid approaching of the bulk limit 7rp(-i)(y) = 1 as y ^ cxd, but slow decay outside the 
boundary of support as 7/ — > — cxd. This observation is easily confirmed analytically from ( ^.33|) . 
Use of the expansion 



1 , e-'^y 
1 + 71^' 



1 + erf (\/2y) 



y -oo 



2n\y\' 



shows 



vrp(i)(y)~ ^ "^'^ (5.35) 

Consider now the density for the OCP with Boltzmann factor proportional to ( |1.7| ). It is 
given by the exact formula ||22| 

^P(i)(r) = e-l^1%(|zp;iV), (5.36) 

so we see from ( ^.320 that in the limit N ^ oo the density as measured from the southern most 
point on the boundary of support is given by 

vr/3(i) (y) = /i(y), h{y) := ^ (l + erf (^/2y)) . (5.37) 



In contrast to the edge density ( 5.34 ) which exhibits algebraic decay as seen in ( 5.35 ) for large 
distances outside the boundary of support, here the asymptotic values are reached with only 
Gaussian corrections in both directions, 



1 , — , , , W — > CXD 



2 



7r/3(i)(y) 



y -oo 



2V2^\y\ ' 



The sharpness of the boundary of support for the density associated with ( [1.7] ), as generated 
from the eigenvalues of complex random matrices, relative to that for the zeros of complex 
random matrices can be illustrated by plotting the zeros for a particular realization in each case. 
This is done in Figure ||. 



Another point of interest is that analogous behaviour to ( |5.35 ) is found for the OCP at 



(3 = 2 with the correct Coulomb potential C + N log r between particles and background outside 
the disk. The exact density as measured from the southern most boundary of the uniform 



background is given by |27] 



(\ ft ^ ^ r exp(23/2tx - 1^) 

P{i){y) = f{y)-j= -r, — f,,. , _t2 r- (5-38) 

^' \/tt Jo 1 + erf(t) + e ' /t0r 
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Figure 1: A typical realization of the eigenvalues of a 81 x 81 complex Gaussian random 
matrix (leftmost plot) and the zeros of a complex Gaussian random polynomial of degree 
81 with variances given by ( |1.6| ). The shaded region represents the disk |z| < 9 which is 
the leading order support of the density in both cases. Outside the disk the density has 
a 1/r^ tail in the case of the zeros, whereas it falls off as a Gaussian for the eigenvalues, 
in keeping with the realizations in the Figure. 



where f{y) = e , y > and f{y) = 1, y < 0. This gives for the asymptotic behaviour [27| 



mi)iy)-{ ' ^W' (5.39) 

An analytic feature of the surface density for the OCP in general is that the total charge is 
zero. This means 

rO 



7rpocp(i)(y) dy+ vr(pocp(i)(y) -'^)dy = 0, (5.40) 

D Jo 

which indeed can be verified for the OCP profiles ( ^.36|) and (|5.3^ ). A straightforward calculation 
using ( p. 33 ) also shows that ( |5.40| ) holds true for the boundary density of the complex zeros. 



5.2 The two-point distribution 



Let /0(2) (yii 2/2; 2;i — X2) denote the N ^ 00 limit of P(2)(n,^2) as specified by ( 3.16 ) with 

zi = xi+i{yi- Vn), Z2 = X2 + i{y2- \^)- (5.41) 
The basic strategy is to explicitly compute the limiting value of the matrix elements of 

U -.= (10 - B^A^B (5.42) 
and the determinant d in ( p. 16 ). From ( |3.15| ) and ( |3.17| ) we have 



d = aiia22 - 012021 = e{\zi\'^;N)e{\z2\'^;N) - \e{ziZ2; N)\'^ 
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With zi and Z2 given by ( ^.41| ) we see that the asymptotics of the first term fohows from ( ^.32| ) , 
while ( ^.311 ) shows 



e{ziZ2]N) ~ ie"i^"2{l + erf(-^(yi+y2-i(3;i-X2))}. 



(5.43) 



Thus, with h{y) defined as in (5.37) we have 



_g-(^i-^2)^-(yi-y2)^ 



(5.44) 



The analysis of the leading order behaviour of the matrix elements of ( [5.42| ) is more involved. 
Each element consists of the sum of many individual elements, and we find that a naive leading 
order expansion of each term individually gives zero, as all such contributions exactly cancel. 
It is therefore necessary to expand individual terms to higher order. This can conveniently be 
achieved by writing the terms e(z; N + p) which occur therein with p 7^ in terms of e(z; A^). 
For example, from the definition (p^) we have 



z^ z^ 
e{z- N-l) = e(z; N) - e{z- N - 2) = e{z; N) - — - _ 

which we use to substitute for e{z; N — 1) and e{z; N — 2). 



(5.45) 



Let us give the details of the implementation of this procedure in the case of uu. First, 
according to ( ^.42 ), we have 



Uu = dcu + 611012621 + 621021611 - 611022611 - 621011621. 
Substituting in the explicit form of the matrix elements according to ( |3.17| ) then gives 

un = (kipe(kiP;^^-2) + e(|zi|2;Ar_i))(e(|zip;Ar)e(|^2|2;Ar)_|e(ziZ2;iV)|') 
-\zi\M\zi\^-,N-l)fe{\z2f;N)-\z2\MziZ2-,N-l)\^e{\zi\^;N) 
+2Re(^ziZ2e{ziZ2; N - l)e{ziZ2; N)e{\zi\'^; N - 1)). 



With the susbstitution ( ^.45 ) this can be rewritten to read 

uu = (a) + (6) + (c) + (d) + (e) + (/) + {gl) + {g2) + (^3), 



where 



(o) := \ei\z^\^■N)\^ei\z2\^■,N) 

(6) := |e(ziZ2;iV)|^e(|zi|2;Ar)(-l- |zi|2-|z2p + ziZ2 + ^2^i) 
(c) := e(|zip;iV)e(|z2p;iV)(-^ + ^(kip-A^)) 



(d) := \eiziZ2;N)\ 



2F1 



Nl 



Nl ^ 

■(\zif + + 1 - Z1Z2 - ZlZ2 
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(e) 
(/) 
(ffl) 
(<?2) 
(53) 



e(|zi|2;iV)e(ziZ2;iV) 



{Z1Z2) 
Nl 
{Z1Z2) 



N 



Z2\ - ZlZ2 



e{\z,\^;N)e{ziZ2;N)^-^^^[\z2\^ - Z1Z2) 



\AN+2 



-e{\z2\^;N) 



|^^|2iV|^2|2JV+2 



(iV!)^ 



-Re(^ 



e(ziZ2;A^) 



Af! V ATI 

It now remains to calculate the large N asymptotic behaviour of each of these terms with zi 
and Z2 given by ( |5.41[ ). Using ( ^.32|) , (|5.43| ) and the expansion 



{ZjZj) 

m 



N 



1 



exp (^N - ViV (yj + yj + i{xj - Xi 



V2ttN 

+ +x'j - yf - + i{x^yi - xjyj] 

„2 / 



a ~ 



which follows from the expansion (1 + u)^ ~ e^^" " /2+--) ^j^^j Stirling's formula, we obtain 

^3N-4VNy.-2VNy, ^^^^2 ^ ^2) + ^2 + yj^^ (y,)h{y2) 

^,N-,VNy.-2VNy, ^ _ _ ^^)2 _ _ ^^^2 _ ^^2 ^ ^2 ^ 2^1X2 + 2yiy2 

x/i(yi) + 2/2 + - 2:2)) 



(d) 

(e) + (/) 



„3Ar-4VlVyi-2VAfj/2 , /I 



yi exp(2xi +xl + yl)h{yi)h{y2) 



~ e 



3iV-4v/F,,-2Viv,, /ly^ exp(x? -yl + 2{x,X2 + yw2)) 
V vr 



^(2^^^ +y2 + «(a:i -2:2)) 



(2(3^1 +yi+A- yl) +xiX2 + yiy2)h{yi) 
2Re|— i=(yi - y2 + iixi - X2))exp {i{x2y2 - xiyi + X2yi - xiy 



. 27r 

x/i(^(yi + y2 - - a;2))) } 



(5I) 
(53) 



3N-.VNy.-2VNy, 1 g^p(2^2 _ ^y2 ^^2^ ^2)^^^^) 



27r 



g3iV-4v^,,-2v^j/2 J_ g^p(2^2 ^ ^2 _ y2)^(y^^ 



27r 



^37V-4v/]V,.-2v^,. 1 g^p (1 (3^2 _ 3^2 + ^2 _ ^2) + ^^^^ + ^^^^ 



27r " V2 

;2Re| exp (i{xiyi - X2y2 + 2:1^2 - 2:2^1 



))/i(^(yi +2/2 - X2)))}. 



Notice that all of the above asymptotic forms contain an N dependent factor e'^^ 4v^?/i 2v^j/2 
Similarly decomposing n22 shows that the asymptotic form of each term contains the depen- 
dent factor e3^-4VlVy2-2VJVs/i_ ^^^^ according to ( |3l5| ) and ( ^ ) 



71" P(2)(n,r2j 



^111^22 , ^12^21 



d3 



+ 



d3 
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so we are required to form the product U11U22 and to divide by d^. Doing this, we see that 
the TV dependent factors cancel out as expected. Similarly we find that N dependent factors 
in ui2U2i/d^ cancel. Collecting together the remaining terms, which are calculated for U22, U12 
and U21 in the same manner as detailed above for uu, we obtain for the limiting value of the 
edge two-point distribution 

,2 . 



, X Vx\V22 , «12'y21 



4 



where 



di := Ky^)h{y2)-e-^-'--^^"-^y'-y^^y{^{yi+y2 + i{x^-X2))) 

- g(3;i-a:2)^+{?/l-2/2)^^_^ 



d2 

Vll 



V22 



V12 



(5.46) 

(5.47) 

(5.48) 



h^{yi)h{y2) - {{xi - X2f + (yi - y2? + l)e 

^(^(yi + 2/2 + i{x\ - X2))) 

'iylp-{xi-x2f-{yi-y2f 



-(x-i-x'i)'^ -{yi-y2f 



Kyi) 



TT 



\{yi + y2 + i{xi-X2))))] (5.49) 



-yie '^y^h{yi)h{y2) 

TT 

/i(^(2/l + y2+iixi - X2 

+2Ref ^(2/1 -y2 + i{xi - x2))e-(f?+f2)e-(^i-^2)V2-(2/i-2/2)V2 
\/27r 

xe'^^^-'=^^^y'+y^^h{yi)h(^^{yi + y2 - i{xi - X2))} 

_ J_e-42/?^(y2) - ^e-^(y'+y'^h{yi) + 2Re|^ 
y^g-3yl-yip-{x2-xi)'^/2-{y2-yif/2^i{xi-x2){yi 

^11 (5-50) 

v\^y2 

-((xi - X2f + (yi - y2f - l)e(^i-^^)'+(^^-^^)'/i(yi)/i(y2)/i(^(yi + 2/2 - ^(xl - xs))) 
-h{^{y\ + 2/2 - «(a;i - ^2))) ^(^(2/1 + 2/2 + ^(a^i - ^2))) 

^(yi + 2/2 + i(xi - X2))e-(2'?+f2)e3(^2-xi)V2+3(2/2-3/i)V2ei(a;i-x2){2/i+2/2)/^(y^);j(y2) 

+ + y2 - ^(ici - x2))e-(2'?+2'i)e(^i-^2)'/2+(2'i-2'2)V2 

^(^(2/1 + 2/2 + i(a;i - 3^2)) 





/2^ 




1 


+- 








xe' 


{xx- 


+- 


1 


\ 


/2^ 


1 


1 





+ 



yi - y2 + ^(xl - X2))e-22/ie(-2--i)'+(^^-^i)'/t(yi)/i(^(yi + y2 - ^(^1 - ^2))) 
y2 - yi + i{x^ - x2))e-2^^e(-^--i)'+(^^-^i)'/i(y2)/i(^(yi + y2 - i{x^ - X2))) 

i_e-2(2/2+j/i)g(x2-xi)2+(2/2-j/i)2g2i(a;i-a:2)(j/i +2/2)^^(^1 + + i{xx - X2))) 
+ ^e-2fe?+2/i)e(-2--i)'+(^2-2'i)'/i(i(yi + 2/2 - i(xi - X2))) 
_J_e-32/f-2/ie3(^2-^i)'/2+3(2/2-2/i)V2e^(^i-=^2)(j/i+yi)/j(y2) 

_ J_e-2/i-32/ie^(^2-a;i)^/2+3(j/2-2/i)V2gi(a;i-j;2)(j/i+2/i);^(^y^^ (5.51) 
27r 
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V21 = Vi2. 



(5.52) 



The above explicit forms show that indeed t^'^ P(2){yi,y2', xi — X2) as specified by (|5.46|) is 
a function of xi — X2 as claimed in the notation. Physically, this corresponds to translation 
invariance in the direction of the boundary. Another check on ( |5.46 ) is that in the limit yi, y2 — > 



cxD, yi — 1/2 fixed, the bulk result ( 3. 181) is reclaimed. In this limit, since erf(t 
we have h{u) — > 1 as u — > cxd. This shows that 



■vii ~ V22 ~ 1 — (r + l)e 



, f 12 ~ V21 ~ 
d2 ~ e"^^ 



1 



1 



1 as M — > 00 



(5.53) 



where := (xi — ^2)^ + (yi — 2/2)^, and hence 



T^'^P{2){yi,y2]Xi - X2) 



(1 



1 + r 



2\2 



which is (|3.1§| ) in the form given by Prosen eq. (27)]. 

It is of interest to compare the structure of ( 5.46| ) with the limiting value of the edge two- 
point distribution function for the OCP with Boltzmann factor ( |1.7D . In the finite system the 
exact two-point distribution is given by [E2I 



IT' 



'p{2)(.ri,f2) = 7rV(l)(n)/9(l)(f2) 



|e(ziZ2;iV)P 



(5.54) 



where 7r/?(x)(r) is specified by ( 5.36 ). Use of (|5.43| ) then gives that in the limit ^ 00 with the 
coordinates (|5.41) we have 



7rV(2)(yi,y2;a;i -X2) = vrV(i)(yi)/5{i)(y2 



■IT e 



-(xi-X2)^-(yi-J/2)^ 



P(i)(^(yi + y2-i(a;i -a;2))) 

(5.55) 

where 7rp(i)(y) is specified by ( |5.37 ). In ( |5.55 ) the general structure of the two-point distribution, 
as a product of the one-body densities plus a correction term (the truncated two-point distribu- 
tion, or two-point correlation function pj^-^) which decays as the points {xi,yi) and (x2,y2) are 
separated, is displayed. Thus we have 

2 



P{2)iyi^y2;xi - X2) 



i)(^(yi 



+ 7/2 - i(Xl - X2 



(5.56) 



This is distinct from the exact expression (5.46) in which such a decomposition is not exhibited. 



5.3 Asymptotics of the two-point correlation 

In this subsection the asymptotic behaviour of the two-point correlation parallel to the boundary 
of support of the density will be investigated. In the case of the two-point correlation ( 5.55| ) for 



the OCP at /3 = 2, this behaviour must conform to a universal form obeyed by the asymptotics of 
the charge-charge correlation parallel to a plane boundary for general two-dimensional Coulomb 
systems in their conducting phase. This asymptotic behaviour in the case of the OCP specifies 
that [^ ] 

PoCP(2)iyi^y2;xi-X2) ~ - -^^^^'^^^ f dyi f ^^2/(^1,^2) = l- (5.57) 

^' |a;i-a::2|-»oo ZpTT^[Xi—X2) J-oo J-oo 
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To test the prediction ( ^.57 ) on the exact result ( 5.55 ), the asymptotic expansion of the error 
function erf(y + iX) for large imaginary argument is required. Now, by choosing the path of 
integration from to y + iX in the definition of erf(y + iX) to first go along the imaginary axis 
to iX then parallel to the real axis to y + iX, we see that it is possible to write 



erf (y + iX) 



e" ds + — 
vr 



Using integration by parts we can then deduce that 

y 



Y 



erf(y + iX) ]^^x--Y--2^XY^ ^ 



\x~Y^ 

This in turn, when used in ( 5.55| ), implies 



PocPi2)iyi^y2;xi - X2) 



/(yi)/(2/2 



fiy) 



d 



(5.58) 



(5.59) 



\xx-x2\^oo 47r2(xi — X2)^ ' '' dy 

which indeed obeys ( |5.57] ). 

Let us now consider the asymptotic expansion of ( 5.46| ) for l^i — X2I 00. Since this 
expression represents the full two-point distribution rather than its truncated counterpart we 
expect the leading behaviour to be given by t^"^ P{i){yi)P(i){y2), where each /0(i)(y) is specified 
by (5.34). The next term will then give the leading order behaviour of p^-^. Now according 
to ( ^.46[) and ( ^.48|) , ( ^.50 ), ( 5.52| ) the asymptotic behaviour follows from the behaviour of the 
quantities di, vu and V12. Inspection of ( 5. 471) , (|5.49|) and ( ^.51 ) shows that the latter asymptotic 
behaviour follows from that of h{\{yi +y2 — i{xi —X2))), its complex conjugate and its modulus 
squared. Now, from the definition ( ^.37| ) and the expansion ( |5.58| ) we have that 

1 



ikyi 



h\ -(2/1 + y2 - i[Xi - X2 



^(^(2/1 + y2 + i{xi - X2) 

+2/2 -i{xi -X2)) 



27r 



62 



( 



yi + 2/2 



+ 



X2-X1 {x2-XiY [X2-X1Y 



r(i-(yi + y2)^) + 



^(^(yi + y2- iixi - X2, 



1 

2^ 

+ 



,{xi-x2f -{yi+y2)'^ 
1 



1 



{Xi - X2) 



V(xi - X2y 
(2-(2/i + y2)') + -- 



After some calculation we then find 

-2(y?+yi) 

d^ 



h{yi)h{y2)[l - „ , . . . . .3 
V 27r/i(yi )h[y2 ) [xi - X2 Y 

h{y2){h\yi) - J^yie~^y"Hy,) - ^e~^y' 



+ 



Vl2 ~ 6 2 



g-2{y2+y2) 

(Xi - X2)^ 
{xi-X2f 



(--%i)(l + 4y?) 



6 2 



2tt V2^ 

{yi+y2)^ ^-i{x2-xi){yi+y2) 



-2yf 



(5.60) 



(5.61) 



X2 - xi ^^/2^^ 



i=4yiy2/i(2/i)/i(y2)e^^^'+^'^ + -yi/i(2/i)e^^i + -y2h{y2)e 

/2-K vr TT 



+ 
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+7 ^(^(yi +y2)(l - 2yiy2)/i(yi)%2)e2(^?+^2') + ^(1 - 2^2(^1 + ^2)) 

y<h{y2)e^yl + ^(1 - 2yi(2/i + y2))h{yi)e^y' - ^(^i + ^2)) } (5.62) 



Substituting the results ( ^.601 ), ( ^.611 ) and ( |5.62| ) in ( |5.46| ), recalling the relations ( ^.48| ), 
(|5.50| ) and ( |5.52|) , and performing some further calculation we find the sought expansion, 



7r^/'(2)(yi,y2;a;i - X2) 



givi) 5(2/2) , Q{yi,y2 



+ 



where 



K^iyi) h'^{y2) {xi - X2Y 
g{y) := h\y) - J'^ye-^v' h{y) - 



(5.63) 



27r 



Qiyi,y2) 



1 



+- 



^((12y2y2 - 3y2)h\yi)hiy2)e^y' + {Uyiy^ - 3yi)hiyi)h\y2)e^y') 



xe-4(ff+s/i) + ^ymHyi)h{y2)e-'''y" 



+ ^((%? - l)/i2(yi)e-2(^'?+3.i) + (4^2 _ i);^2(y^)^-2(3.H.i) 

+ ^(6yi/j(yi)e-2(2,?+3yi) ^ 6y2/i(y2)e-'(='^?+'^2') 



_L^p-6(y?+yi) 



(5.64) 



Comparison with ( 5.34 ) shows that indeed the leading order behaviour is t^'^ P{i){yi)P{i){y2) so 
we have 

Pi2M,y2;x,-x2)r^^^^^^^^--^. 

Furthermore, a close inspection of ( |5.64| ) and ( ^.34] ) reveals that it factorizes as 

Q{yi,y2) = f{yi)f{y2), f{y) ■= ^^/'(i)(y)- 



Consequently we have 



P(2)iyi^y2;xi - X2) 



f{yi)f{y2) 



f{y)dy = i, 



(5.65) 



47r2(xi — X2y ' 

which is identical to the behaviour ( |5.57| ) for the OCP at /3 = 2, except for the sign. 

Note also that since for small separation between {xi,yi) and (a;2,y2) the truncated distri- 
bution is negative, the positive tail for large |rEi — X2I exhibited by ( 5.57| ) indicates that pj^^ 
changes sign, analogous to its behaviour in the bulk. In contrast the same quantity for the OCP 
at /? = 2 is always negative, as seen from ( ^.56 ). 
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6 Summary 



The p.d.f. (|l.3D when interpreted as a Boltzmann factor consists of the pairwise logarithmic 
interaction (1.5) as well as an extensive many body interaction. If the extensive many body in- 
teraction is replaced by a one body Gaussian factor, the p.d.f. becomes identical to the p.d.f. (1.7) 
giving the eigenvalue distribution for complex Gaussian random matrices, or equivalently the 



Boltzmann factor for a 2dOCP. The exact distribution functions for both ( |1.3| ) and (1.7) can be 
calculated exactly, enabling a comparative study of the two systems to be undertaken. 

Some similarities, and some differences, have been found. The similarities include the same 
y ^ —oo decay of the density profile away from the boundary being displayed by ( |5.35| ) for 
the complex zeros, and by ( 5.39D for the soft wall OOP. Also both systems have the same (up 
to a minus sign) universal form ( |5.57| ) of the two-point correlation parallel to the boundary. A 
noteable difference is that while the second moment of the two-point function in the bulk of 



the OOP obeys the Stillinger-Lovett sum rule (4.9), the same quantity for the complex zeros 
vanishes. This was shown to have consequence regarding the behaviour of the variance of a 
linear statistic, a quantity which has received recent attention in the context of two-dimensional 
Coulomb systems |17], and the zeros of multidimensional random polynomials with a constant 
density on the surface of higher dimensional spheres |14]. 
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